function xdeltapcom_sca = NAO_deltapcom_sca(x,  p, a)
xdeltapcom_sca = (-887438959316366215603613875869840 - 213282321474164819037149660846989289* ...
   x(1) - 125924008303117708982302087981441771*x(2) -  ...
  47361140992860211861416779584211771*x(3) -  ...
  25518317987395192388327226340538229*x(4) -  ...
  9655073790944689329642589190113229*x(5))/ ...
 1061120650506889946073703089833150000 ...
;